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The fluxes of energetic particles in the radiation belts are found to be strongly controlled by the solar wind 
conditions. In order to understand and predict the radiation particle intensities, we have developed a 
physics-based Radiation Belt Environment (RBE) model that considers the influences from the solar wind, 
ring current and plasmasphere. Recently, an improved calculation of wave-particle interactions has been 
incorporated. In particular, the model now includes cross diffusion in energy and pitch-angle. We find that 
the exclusion of cross diffusion could cause significant overestimation of electron flux enhancement 
during storm recovery. The RBE model is also connected to MHD fields so that the response of the radiation 
belts to fast variations in the global magnetosphere can be studied. We are able to reproduce the rapid flux 
increase during a substorm dipolarization on 4 September 2008. The timing is much shorter than the time 
scale of wave associated acceleration. 

Published by Elsevier Ltd. 


1. Introduction 

The Earth’s radiation belts consist of energetic electron ( ~ 100 keV 
to several MeV) and ions ( ~ 1 00 keV to several hundred MeV) trapped 
in the magnetosphere roughly from 1 .2 < L < 8. The energetic elec- 
trons reside in 2 distinct regions: the inner belt and the outer belt, 
which are usually separated by the slot region ( 1 .8 < L < 3 ) of depleted 
particle populations. Pitch-angle diffusion loss of electrons by inter- 
acting with whistler mode plasmaspheric hiss is believed to be the 
cause of the slot region (Lyons et al„ 1 972 ; Albert, 1 994; Meredith et al„ 
2007). The inner belt is relatively stable while the outer belt is highly 
variable with geomagnetic activity. The fluxes of energetic electrons in 
the outer belt decrease during the main phase of a magnetic storm due 
to adiabatic effect (Dessler and Karplus, 1961; Kim and Chan, 1997). 
Additional non-adiabatic processes also contribute to the flux decrease 
in the storm main phase (Green et al., 2004; Ukhorskiy et al., 2006). 
During the recovery phase the flux of energetic electrons can change 
dramatically as well. While approximately half of all moderate and 
intense storms cause a net increase in the flux of energetic electrons by 
a factor of 2 or more, approximately a quarter of these storms result in 
a net decrease in the fluxes by more than a factor of 2 (Reeves et al„ 
2003 ). This variability is caused by an imbalance between acceleration, 
transport, and loss processes all of which become enhanced during 
geomagnetic storms (Horne, 2002; Thorne et al., 2005; Horne et al„ 
2006; Summers et al., 2007). The ratio of post-storm to pre-storm flux 
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was found related to the solar wind speed (Paulikas and Blake, 1979; 
Reeves et al., 2003) and the direction of the IMF Bz during the storm 
recovery phase (lies et al., 2002). On the other hand, Ukhorskiy and 
Sitnov (2008) suggested that the outer belt can respond differently to 
similar solar wind driving. There is a broad range of processes that 
shape the radiation belts. Some of them are nonlinear mechanisms of 
local particle acceleration, such as, interactions with whistler mode 
chorus waves (Summers et al., 2004; Omura and Summers, 2006), and 
relativistic electrons drift-resonance with large amplitude, narrow 
bandwidth Pc 5 waves (Degeling et al„ 2008). 

The intensification of the radiation belts has significant space 
weather consequences. Moderate energy ( ~ 1 0- 1 00 keV) electrons 
can cause surface charging effects and relativistic (~0.1-5 MeV) 
electrons can cause deep-dielectric charging on space systems 
(Baker, 2001 ). Therefore, understanding the physical processes that 
are controlling the development of the radiation belts during active 
periods and being able to predict their variability have both 
scientific and practical significance. 

There exist multiple sources of radiation belt particles. Radial 
diffusion has traditionally been considered to be the leading 
transport and energization mechanism in the inner magnetosphere 
(Schulz and Lanzerotti, 1974). However, it has recently been 
suggested that electrons can be accelerated efficiently by 
resonant wave-particle interactions with whistler mode chorus 
waves (Horne and Thorne, 1998; Summers et al., 1998) and fast 
magnetosonic waves (Horne et al., 2007). 

A number of kinetic models have been established to simulate 
the radiation belt dynamics and to provide interpretation for 
observable features. In a kinetic model, the equation for the particle 
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distribution function is solved analytically or numerically. One 
simple approach is based on a standard radial diffusion equation 
with diffusion coefficients driven by the solar wind conditions or 
geomagnetic activity (Li et al., 2001; Albert et al., 2001). Varotsou 
et al. (2005) and Horne et al. (2006) combined radial diffusion with 
acceleration and loss due to whistler mode chorus waves and 
confirmed that wave acceleration by whistler mode chorus is an 
important acceleration mechanism in the outer radiation belt. By a 
similar approach but including cross diffusion terms due to chorus 
waves, Albert et al. (2009) simulated the 9 October 1990 magnetic 
storm. They found both chorus wave acceleration and radial 
diffusion were required to account for the observed increase in 
relativistic electron flux during the recovery phase. For models that 
cover a wide range of energy, drift motion must be considered since 
convection is an important transport mechanism for lower-energy 
( < 50 keV) particles (Bourdarie et al„ 1997; Zheng et al., 2003; 
Miyoshi et al., 2006; Fok et al., 2008). Using the relativistic kinetic 
model with dipole magnetic field, Miyoshi et al. (2006) simulated 
the dynamics of energetic electrons during the October 2001 storm. 
They reproduced the observed local time flux asymmetry for hot 
electrons (30 keV). They also concluded that only convective 
transport and radial diffusion cannot explain the enhancement 
of relativistic electrons seen during storm recovery. 

A convection-diffusion model, namely the Radiation Belt Envir- 
onment (RBE) model, has been developed to understand the 
radiation belt dynamics in order to predict the flux variation 
during active times (Fok et al., 2001, 2005, 2008; Zheng et al., 
2003). The RBE model employs time-varying, realistic magnetic 
field so that radial diffusion effects due to slow magnetic 
fluctuations and inductive radial transport owing to fast 
magnetic reconfiguration can be modeled. The RBE model was 
used to simulate a substorm injection during a dipolarization of the 
magnetic field (Fok et al., 2001). Observable features during 
substorms, such as dispersionless injection and drift echoes, are 
successfully reproduced. Electron flux enhancements during 
magnetic storms were also studied using the RBE model (Zheng 
et al., 2003; Fok et al., 2005, 2008). They found that energization by 
the inductive electric field and by whistler mode waves is crucial 
for the flux increase during magnetic storms. A simplified version of 
the RBE model is currently running in real-time to provide radiation 
belt now-casting updated every 15 min. The geosynchronous 
fluxes at longitudes of GOES-11 and 13 are extracted from the 
RBE real-time run and are plotted together with real-time GOES 
electron ( > 0.6 MeV) data. The model-data comparison is 
continually posted at http://mcf.gsfc.nasa.gov/RB_nowcast/. 

Recent developments, discussed in this paper, expand on past 
work in a number of ways. In previous RBE simulations of wave- 
particle interactions, only pure energy and pure pitch-diffusion are 
considered (Fok et al., 2005, 2008). Since cross diffusion may have a 
comparable effect to pure diffusions (Albert and Young, 2005), the 
RBE model has been extended to include cross diffusion in energy 
and pitch-angle. We choose to use the Alternating Direction 
Implicit (ADI) method (Xiao et al., 2009) to solve the cross 
diffusion terms in the RBE equation. We revisit the geomagnetic 
storm on 23-27 October 2002 and quantitatively assess the 
influence of cross diffusion in electron flux enhancement during 
the recovery phase. Another development in the RBE model is 
connecting the model with a global magnetohydrodynamic (MHD) 
model. Glocer et al. (2009) coupled the RBE model with the Block- 
Adaptive-Tree-Solarwind-Roe-Upwind-Scheme (BATSRUS) MHD 
model to simulate radiation belt development during steady 
solar wind conditions and during a real storm event. They found 
the coupled code was able to create rapid inward transport on the 
time scale of tens of minutes. In this work, we use the RBE-BATSRUS 
code to simulate a substorm injection during the 3-5 September 
2008 storm. We examine the timing and location of electron 


enhancement during the dipolarization event. In the following, a 
brief description of the RBE model is given. We then outline the 
algorithm of the ADI method and present the simulation results 
with cross diffusion included. Finally we present the RBE simula- 
tions of a substorm dipolarization and compare the calculations 
with Akebono electron data. 


2. Radiation belt environment (RBE) model 


The RBE model is a kinetic model that solves the bounce- 
averaged Boltzmann equation to obtain the distribution function of 
energetic electrons (Fok et al., 2008) 
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where 

G = T(a 0 )sin2a 0 (E-|-£ 0 )v / £(£+2£ 0 ) 

/= ( t, X,, </>;, M, K ), is the average distribution function on the field 
line between mirror points. X t and </>,- are the magnetic latitude and 
local time, respectively, at the ionospheric foot point of the 
geomagnetic field line. M is the relativistic magnetic moment 
and K = J/J&m 0 M, where J is the second adiabatic invariant. The 
motion of the particles is described by their drifts across field lines, 
which are labeled by their ionospheric foot points. The inner 
boundary of /,■ is at 11.8 , corresponding to L=1.06. The outer 
edge of the model is bounded by field lines with X, not greater than 
70.2 ' and an equatorial crossing at 10 Earth radius (R £ ), whichever 
is closer. The M range is chosen to fully represent the energies of 
electrons from lOkeV to 6 MeV. The K range is chosen to cover 
the loss cone so that particle precipitations can be estimated as 
well. £ 0 is the electron rest energy, a 0 is the equatorial pitch-angle, 
and T(a 0 ) is given by 

T(a 0 )= 1 T'"— (2) 

R 0 Jo cos a 

where R a is the radial distance in R E of the field line equatorial 
crossing. The integration is along the field line from the equator to 
the mirror point. x b is the particle bounce period. 

The left hand side of (1) represents the drifts of the particle 
population and the terms on the right hand side of (1) refer to 
diffusion and loss. The calculation of the bounce-averaged drift 
velocities across field lines, <2,> and <</>,), were described in 
detail in Fok and Moore (1997). These drifts include gradient and 
curvature drift, ExB drift from convection, and corotation electric 
fields. In this "drift" approach, the RBE model can simulate the 
effects of radial diffusion only if the electric and magnetic field 
fluctuations are properly and consistently represented. 

The effects of the inductive electric field due to a time-varying 
magnetic field are also taken into account implicitly in the model 
(Fok et al., 2005). We have assumed that field lines are rooted at the 
ionosphere, so that the inductive electric field there is zero. 
However, the shapes of field lines at higher altitudes vary as a 
function of time according to the magnetic field model. If field lines 
are perfect conductors, the field line motion at high altitudes, e.g., 
at the equator, will generate an induction electric field (E md ) of 
the form, 

Eind = — Vo X B 0 (3) 
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Fig. 1. An illustration to show the generation of the inductive electric field (E inli ) 
during field-line re-configuration. 

where v D and B 0 are the field line velocity and magnetic field at the 
equator. Fig. 1 is an illustration on field line mapping changing with 
time and the generation of E ind . Under the frozen-in condition, 
particles initially moving along a particular field line at G will 
continue to share the same field line at t 2 when the field line 
becomes more stretched, as shown in Fig. 1 . Particles are drifting 
outward and decelerated in this case, and the first and second 
adiabatic invariants are conserved if the field line motion is slow. In 
the RBE model, electrons are energized or de-energized during the 
course of varying the magnetic configuration. The electron pitch 
angles will change as well. 

The first two terms on the right hand side of (1) represent 
particle diffusion in energy and pitch-angle from interactions with 
plasma waves. When solving these two terms, we first map the 
particle phase space density from to ( E,a 0 ) coordinates, 

perform diffusion in E and a 0 , and then map the updated 
distribution back to the ( M,I < ) coordinates (Fok et al., 1996). The 
diffusion terms are followed by losses due to the loss cone, the 
boundary of which is assumed to correspond to mirror height of 
120 km. Particles in the loss cone are assumed to have a lifetime of 
one half bounce period (0.5 r b ) (Lyons, 1973). 

Eq. (1) includes multiple processes of different timescales. 
We use the method of fractional step or operator splitting to 
decompose the equation and solve only one term at a fractional 
step (Fok et al„ 1993). To solve (1), we have to specify the electric 
and magnetic fields, the initial distribution, and the particle 
distribution on the nightside boundary, which is set at 10 R E or 
the last closed field line. The effect of radial diffusion is 
incorporated via these time-varying electric and magnetic fields. 
The NASA trapped radiation model (AE8MAX) (Vette, 1991; Fung, 
1996) is used for the initial condition in the entire RBE spatial 
domain. The distribution at the nightside boundary is assumed to 
be a kappa function with density ( N ps ) and characteristic energy 
(E ps ) modeled by linear relations with the upstream solar wind 
conditions (Zheng et al., 2003) 

Nps(t) = [0.02N SW ( t-2 h) + 0.3 1 6] v/amu 

Eps(t) = 0.016V sw (t— 2h)— 2.4 (4) 

where N ps is in cm~ 3 , N w is the solar wind density in the same unit, 
amu is the atomic mass unit of the electron, E ps is in keV, and is 
the solar wind velocity in km/s. Note that we assume a 2 h time lag 
between the plasma sheet condition and solar wind condition at 
the dayside magnetopause (Borovsky et al., 1998). 

3. Cross diffusion in the RBE model 

We have previously simulated the evolution of radiation belt 
electrons during the 23-27 October 2002 geomagnetic storm using 
the RBE model (Fok et al„ 2008). The magnetic field is specified by 
the Tsyganenko 2004 (T04) model (Tsyganenko et al., 2003) and 


electric field by Weimer 2000 model ( Weimer, 2001 ). Magnetic and 
electric fields are updated every 5 min according to the 
instantaneous solar wind condition and Dst index. The bounce- 
averaged diffusion coefficients are given by the Pitch-Angle and 
Energy Diffusion of Ions and Electrons (PADIE) code (Glauert and 
Horne, 2005). Only resonance with lower-band whistler mode 
chorus (0.1 f ce <f< 0.5 f ce ) is considered. The presence of chorus 
waves is confined between-15° and 15° magnetic latitude. The 
diffusion coefficients are calculated as a function of L shell, energy, 
pitch-angle and / pe // ce , the ratio of plasma frequency to the 
cyclotron frequency. To calculate / pe //ce. which depends on the 
plasma density (n e ), we have embedded inside the RBE model 
the core plasma model of Ober and Gallagher (Ober et al., 1997). 
The PADIE diffusion coefficients are scaled with a chorus wave 
intensity of 10 4 pT 2 . To obtain the actual diffusion coefficients, we 
estimate the chorus intensity at a given location and time during 
the storm using the survey of CRRES plasma wave data for lower- 
band chorus presented by Meredith et al. (2001, 2003). For our 
application the wave data were binned in L shell, magnetic local 
time and 3 levels of magnetic activity (Kp < 2; 2 < Kp < 4; Kp > 4). 

In Fok et al. (2008), we only considered pure energy and pure 
pitch-angle diffusion. Recently, we have implemented cross 
diffusion in the RBE model. The cross diffusion terms can be 
solved using various methods or numerical schemes (Albert and 
Young, 2005; Tao et al., 2008; Xiao et al., 2009; Shprits et al., 2009). 
We choose the Alternating Direction Implicit (ADI) method, 
which we found is stable, efficient and easy to implement. In the 
ADI scheme, the cross diffusion terms in (1) are rewritten as (Xiao 
et al., 2009) 

!-ei( GD -*KA( CD '-£) 

= 1 8(GD) 8f f 1 d(GD) 8f ^ ^ d 2 / 

G 8a. 0 8e C d£ 8a. 0 da. 0 d£ 

=hf+y+hf (5) 

where 

e = E/E 0 , D = eDexJE = eD^e/E 

The cross diffusion terms are thus expressed as derivative in e 
(fi), derivative in a D (f 2 ) plus mix derivative in e and a 0 (f 3 ). To 
advance /in time, a time step is divided into two sub-steps as 


^[hf n+y2 +i2f n +hr] 

(6a) 

■+^[hf n+V2 +y n+ ' > +hf n+v2 ] 

(6b) 


where 

f n =f(t) att = nAt 

In the first sub-step, G is solved implicitly and I 2 explicitly. In 
turn, in the second sub-step, I 2 is solved implicit and G explicitly. f 3 
is always solved explicitly. The finite difference representation of 
(6) is given in the Appendix. 

Having previously simulated the energetic electron fluxes during 
the 23-27 October 2002 event (Fok et al., 2008), we revisit this event 
with cross diffusion terms included in the governing equation. Fig. 2 
shows the L-time diagrams of equatorial electron fluxes in energy 
ranges of 20-70 keV (left panels) and 0.6-1. 8 MeV (right panels); 
radial transport due to time-varying magnetic and electric fields are 
included in all three rows. The top panels of Fig. 2 display results in 
which wave diffusion is not included in the calculations. The Dst 
index (black curve) is overlaid on the plots. The ring current electrons 
(Fig. 2a) and MeV electrons (Fig. 2b) behave very differently in the 
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23 24 25 26 27 23 24 25 26 27 

October 2002 October 2002 

Fig. 2. Simulated electron fluxes on 23-27 October 2002. Left panels: 20-70 keV. Right panels: 0.6-1 8 MeV. Top panels are fluxes without wave-particle interactions. Middle 
panels are fluxes with energy and pitch-angle diffusion. Bottom panels are fluxes with cross diffusion included. The black curves in the top panels are Dst. 


storm. During the main phase, ring current electrons from higher L 
shells drift earthward and fill the entire outer belt. The fluxes remain 
high because no wave associated loss mechanism is included in this 
calculation. For MeV electrons (Fig. 2b), magnetic field effects 
dominate over convection. Noticeable flux dropout is seen in the 
heart of the outer belt during the main phase when ring current is 
intensified. This is the well known Dst effect (Dessler and Karplus, 
1961 ; Kim and Chan, 1997). The magnetic field produced by the ring 
current inflates the main field. The drift shells of energetic particles in 
the radiation belt expand outward correspondingly (conservation of 
the third adiabatic invariant) and particles decelerate. This de- 
energization causes an adiabatic decrease in particle fluxes since 
the energy spectrum slope is negative in the radiation belt energy 
range. The expansion of drift shells also produces permanent loss at 
high L’s when particles encounter the magnetopause (Ukhorskiy 
et al., 2006). During the recovery phase, MeV electron fluxes 
recover as well. The intensities are higher than the pre-storm level. 


Fok et al. (2008) found that the flux increase is a result of electron 
injection and earthward transport during the course of the storm. 

Next we examine the effects of wave-particle interactions on 
the ring current and radiation belt electrons. The middle panels of 
Fig. 2 show the RBE electron fluxes with pure energy and pitch- 
angle diffusion from interacting with whistler mode chorus. 
Calculation results with cross diffusion are plotted in the bottom 
panels. At the beginning of the storm, wave activity is weak. There 
is no obvious difference between runs with and without chorus 
waves. During the main and early recovery phases, wave activity is 
strong. For ring current electrons (left panels), the effects of pitch- 
angle diffusion are stronger than those of energy diffusion. 
A significant amount of electrons are diffused into the loss cone. 
A flux hole is formed around L=4 during the storm recovery. In 
contrast, with the consideration of chorus waves, MeV electron 
fluxes gradually increase during the recovery phase (right panels, 
Fig. 2). The injection of lower-energy electrons forms a seed 
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population for energy diffusion. Particles are accelerated by the 
chorus waves and MeV electrons slowly diffuse to larger L shells at 
late recovery. 

When comparing the fluxes in middle panels with those in the 
bottom panels of Fig. 2, we find noticeable effects of cross diffusion. 


flllX (D EH+Dooao) / flux(DKK+Daoao+DKao) 



Fig. 3. Flux ratio as a function of time of MeV electrons without and with cross 
diffusion. Ratios are calculated at simulation time of 4 days. 


Cross diffusion tends to moderate the impacts from pure pitch- 
angle and energy diffusion. For ring current electrons, with cross 
diffusion included (Fig. 2e), pitch-angle diffusion is weakened and 
the overall fluxes are higher than those without cross diffusion 
(Fig. 2c). Similarly in the case of MeV electrons, cross diffusion 
reduces the energy gain by energy diffusion and thus the flux 
enhancement in the recovery phase. Fig. 3 illustrates quantitatively 
the effect of cross diffusion in radiation belt electrons. The ratio of 
MeV electron flux without and with cross diffusion is plotted as a 
function of L shell at the end of 4-days simulation. In most parts of 
the inner magnetosphere, the ratio is close to 1. However, in the 
heart of the outer belt around L— 4, ignoring cross diffusion could 
cause overestimation of electron flux as much as a factor of 5. 
Detailed analysis of the effects of cross diffusion will be reported in 
a separate study. 


4. RBE simulation of a MHD substorm 

The Sun was in its deep minimum in years 2008 and 2009; there 
was minimal solar and geomagnetic activity. However, recurring 
high speed streams could trigger substorms in the Earth’s magneto- 
sphere even in solar minimum (Baker et al., 1998). On 3 September 



Fig. 4. Dst, AU, AL and solar wind speed, density, By and Bz on 3-6 September 2008. 
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3 Akebono Electrons: >2.5 MeV 



September 2008 

Fig. 5. Top panel: L-Time plot of Akebono electron flux ( > 2.5 MeV) on 3-4 
September 2008. The black curve is Dst. Middle panel: corresponding RBE simulated 
flux with T04 magnetic field. Bottom panel: simulated flux with BATSRUS fields. 


2008, a high speed stream arrived at the magnetosphere and 
triggered a moderate storm with minimum Dst of -51 nT on 
September 4, 0500 UT. Fig. 4 plots the Dst, AU, AL and solar wind 
speed, density, By and Bz on 3-5 September 2008. 

The Akebono satellite observed enhancements of radiation belt 
electrons during the active period on 3-5 September 2008. 
Akebono was launched in February 1989 by the Institute of Space 
and Astronautical Science in Japan (Takagi et al., 1993). In 
September 2008, Akebono was in a high inclination, highly 
elliptical orbit with apogee at 5260 km altitude, perigee altitude 
at 295 km, and orbit period of 2.5 h. The Radiation Monitor (RDM) 
measured electron fluxes in three energy channels : 0.30-0.95 MeV, 
0.95-2.5 MeV and > 2.5 MeV (Takagi et al., 1993). A rapid increase 
in high-energy electron flux with a factor of 80 at 4 < L < 5 is seen 
between 0321 and 0546 UT, around the peak of the storm 
(minimum Dst). Similar growths are also seen in the other 2 
RDM lower-energy channels. The top panel of Fig. 5 shows the 
Akebono electron flux of energy > 2.5 MeV from September 3, 
00 UT to September 4, 12 UT. The Dst index during this time is 
overlaid in the plot. The Akebono data are averaged over 3 orbit 
periods. The time scale of this enhancement is a few hours and thus 
it is too short for wave associated energization. The large AL values 
in this time period (Fig. 4) suggest a substorm may play a role in this 
quick increase of outer belt electrons (Nagai et al., 2006). 

We perform RBE simulations to understand the cause of the flux 
enhancement seen in the Akebono electron data. In most of our 
previous RBE calculations, we used empirical models of magnetic 
and electric fields (Fok et al., 2001, 2008; Zheng et al., 2003). 
Recently, we have adapted the magnetic and electric fields output 
from the coupled model of BATSRUS-Rice Convection Model (RCM) 
into the RBE model in order to self-consistently simulate the 
responses of the radiation belts to solar wind and ring current 
variations during storm time (Glocer et al., 2009). The MHD 
electromagnetic fields are updated every 10 s to drive the drift 
motion and radial transport of radiation belt particles. Fig. 5(b) 
depicts the RBE electron fluxes calculated with T04 magnetic field 
model and Weimer electric field model (Tsyganenko et al., 2003; 
Weirner, 2001 ). Fig. 5(c) is the RBE flux calculated in the MHD fields 
simulated from the BATSRUS-RCM model (De Zeeuw et al., 2004). 
Note that the RBE fluxes shown in Figs. 5(b) and (c) are equatorial 
fluxes and with temporal resolution of 1 h. In contrast, the Akebono 
measurements are taken along high inclination orbits. However, 
the temporal variability of high latitude fluxes was found nearly 
identical with the equatorial fluxes (Kanekal et al., 2001, 2005). 
Wave-particle interactions are not included in these RBE 
calculations. In the quiet period on September 3, the two RBE 
simulations give similar flux intensity. During the main phase of 
the storm, both RBE runs produce flux dropout in the outer belt, in 



Fig. 6. BATSRUS-RCM simulation at (a) 05:00 UT and (b) 05:20 UT on 4 September 2008. Magnetic field lines (white lines) and pressure (color) are plotted on X-Z plane. 
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response to the ring current intensification simulated in the T04 
and BATSRUS-RCM models. However, at 05-06 UT on September 4, 
a sudden increase in electron flux is seen at 4 < L < 5 in the RBE- 
BATSRUS calculation, consistent with the Akebono data. There is no 
significant enhancement in the RBE-T04 run during the recovery 
phase of the storm. 

Since wave-particle interactions are not considered in these 
particular RBE calculations, the enhancement seen in Fig. 5(c) must 
be a result of particle transport. When we examine the magnetic 
configuration during the enhancement, we find the MHD model 
predicts a substorm dipolarization at ~05 UT on September 4. Fig. 6 
shows the BATSRUS field lines in white and pressure in color on the 
X-Z plane before (left) and after (right) the substorm onset on 
September 4. At 05 : 00 UT, the field lines on the nightside are greatly 
stretched. 20 min later, the dipolarization takes place in the tail. 
Field lines are converting earthward and have more dipole-like 
shape. Electrons which are gyrating along these collapsing field 
lines can be accelerated significantly on a time scale of minutes (Fok 
et al., 2001 ; Glocer et al., 2009), much faster than the time scale for 
energization by whistler mode chorus waves, which is typically of 
the order of 1-2 days (Summers and Ma, 2000; Horne et al., 2005). 
The T04 model, which is driven by Dst and solar wind parameters, 
does not contain clear substorm signatures. Empirical magnetic 
field models of this kind cannot directly simulate substorm 
reconfiguration unless special tricks are applied (Delcourt et al., 
1990, 1997; Pulkkinen et al„ 1991; Fok et al., 2001). For moderate 
storms such as this one on 3-5 September 2008, convection is weak 
and the dominant energization and transport mechanism is sub- 
storm reconfiguration and the resulting dipolarization electric field. 
The rapid enhancements of radiation belt fluxes during modest 
storms cannot be explained without the consideration of substorm 
effects. 


5. Discussion and conclusions 

One must be cautious to interpret the timing signature seen in 
the Akebono data. The temporal resolution is limited by the orbit 
periods (2.5 h). In our separate study of the storm in September 
2008, we will identify storm and substorm signatures using 
continuous high resolution data, such as measurements from 
NOAA and GOES satellites. 

So far we have considered only low-latitude whistler mode 
chorus waves. However, study has shown that chorus waves at high 
latitude are important in the loss and acceleration of energetic 
electrons (Horne and Thorne, 2003). Other wave modes, such as 
plasmaspheric hiss and electromagnetic ion cyclotron waves, also 
play crucial roles in the development and decay of the radiation 
belts (Meredith et al., 2007; Lorentzen et al., 2000; Summers and 
Thorne, 2003). All these wave modes have different sources of 
excitation and are found in different regions of the inner 
magnetosphere. We plan to gradually include all the important 
wave modes in the RBE model. In that case we will be in a better 
position to understand and identify the physical processes that 
control the observed variability in the radiation belts. 

Inclusion of the cross diffusion term does make a noticeable 
difference in electron flux in the heart of the outer belt. However, 
the effect is relatively mild and localized when compared with 
other processes such as particle injection, transport and accelera- 
tion. It is difficult to show quantitatively that with cross diffusion 
will improve the data model comparison. Nevertheless, inclusion of 
cross diffusion will give a better estimation of the diffusive effect 
from interacting with a particular wave mode. 

The dynamics of radiation belt electrons is strongly controlled 
by the magnetic configuration and its fluctuations. The empirical 
models of Tsygenenko (Tsyganenko and Stern, 1996; Tsyganenko 


et al., 2003) have been used to simulate the temporal variations of 
magnetic field in the RBE model. We found strong electron 
energization by the inductive electric field associated with the 
time-varying magnetic field (Zheng et al., 2003; Fok et al., 2008). 
However, as shown in Section 4, we cannot reproduce the rapid flux 
enhancement during a substorm dipolarization with the T04 
model. In contrast, the RBE model driven by MHD fields 
successfully produces the observed increase in electron flux 
during a substorm. Another limitation of the combined RBE-T04- 
Weimer models is that they cannot fully simulate the effect of 
radial diffusion even with T04 and Weimer model input parameters 
updated in time. Both the Tsyganenko and Weimer models are 
empirical and not consistent with each other. This problem can be 
resolved by, again, applying MHD fields in the RBE model. Recently 
work by Huang et al. (2010) has shown that the ULF waves (mHz 
range) predicted by the Lyon-Fedder-Mobarry MHD model well 
represent the ULF wave data observed by the GOES satellites. 

In summary, we have reported recent developments and 
improvements in our radiation belt model. We have implemented 
pitch-angle-energy cross diffusion in our wave diffusion calcula- 
tion. We have simulated a substorm dipolarization event with the 
RBE model embedded in MHD magnetic and electric fields. The 
findings from this model development work include 

(1) Cross diffusion moderates the effects of pure pitch-angle and 
pure energy diffusion. Exclusion of cross diffusion would 
significantly overestimate the flux enhancements of relativistic 
electrons during storm recovery. In our simulation of the storm 
on 23-27 October 2002, at the heart of the outer belt, the 
overestimation can be as high as a factor of 5. 

(2) The strong inductive field during substorm dipolarization 
produces rapid increase in energetic electron flux on a time 
scale of an hour, much shorter than that from wave 
acceleration. 


Acknowledgment 

This research was supported by NASA Science Mission Direc- 
torate, Heliophysics Division, Living With a Star Targeted Research 
and Technology Program, under Work Breakdown Structures: 
936723.02.01.06.78 and 936723.02.01.01.27. 


Appendix: Finite difference representation of the ADI scheme 


Eq. (6) outlines the ADI scheme in solving the cross diffusion 
equation. The second order finite difference discretization of (6) can 
be written as 


1 f n + l/2 f n + l/2 i f-n + 1/2 

'rt'k.m F l k,m>k 


Jk + 1 ,m 


— fk,m + ^k,m(fk,m + 1 ~fk,m 

+ /4,m (fk+ l,m+ 1 + fk - 1 ,m- 1 ~fk - 1 ,m + 1 ~fk+ 1 ,m - 1 ) 


(7a) 


,,2 rn + 1 ,fi+l_,,2 rn + 1 

_ f n + l/2 i { f n+t/2_ f n+t/2\ 

~Jk,m “ r F l k,in \Jk + l,m J k-t ,m ) 

4 _,,3 ( f n + 1/2 f n + 1/2 _ f n + l/2 _ f n + l/2 \ 

~ t F l k,m (Jfc+l.m + l ~<~Jk-t,m-l Jfc-l.m + 1 Jk+t.m-tJ 

where 

i _ At Cfc m+1 D) iim+1 ~Gj fim _ 1 Dfc im _i 

2 Gk m (ftom + l ~&om 1 )(£fc + l — ) 


(7b) 
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2 _ At Gk+-[,mDk+'l,m—Gk-\,mDk-'l,m 

2 Gfc m (Ot om + i — OC om _i )(£j( + 1 

3 = AtD k>m 

(^om + l — ^om-1 )(%+ 1 — £fc-l) 

k is e index and m is a 0 index. 

Eq. (7) represents two tri-diagonal systems similar to the 
Crank-Nicolson method. There are well established numerical 
techniques and stability analysis for this type of problem 
(LeVeque, 2002; Burden and Faires, 2004). 
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